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Abstract. Atmospheric muon neutrino in Deep Core (whose rate and spectra might be soon 
available) should exhibit an ^^anomaly^^ (due to tens GeV up-going muon neutrino suppression, 
converted into tau flavor) that must be imprinted in out-coming rate spectra. We estimate here 
our independent muon neutrino spectra based on SK and its projected record on Deep Core 
Channels. Our estimate (based on cosmic rays, muon records and tested Super-Kamiokande 
(SK) data) differs both in shape and in rate from other previous published spectra. The expected 
rate might exhibit a minimum near channel 6 of Deep Core strings and it should manifest strong 
signature for flavor mixing (mostly between channel 4-15) and a relevant ^^anomaly^^ for eventual 
CPT violation (MINOS like) written at channel 3-6, whose statistical weight (mainly at channel 
5) might soon confirm or dismiss MINOS CPT claim |2j. At the flux minimum around channel 
6, (a flux suppressed respect the non oscillated case at least by an order of magnitude) the 
atmospheric neutrino paucity offers a better windows to a twenty GeV Neutrino Astronomy. 
Therefore by doubling the string array we may foresee a richer rate and a more complete (zenith 
and azimuth) atmospheric neutrino distribution and an exciting first twenty GeV Astronomy. 



1. Introduction 

The role of neutrino muon and tau mixing among flavors has been observed since nearly 13 years 
and it has been foreseen already in different recent high energy astrophysical roles [T],fn],[TU]. 
Recently MINOS results hint for a CPT violation in neutrino sector. We considered therefore a 
long base experiment to Deep Core to disentangle such a neutrino CPT puzzle. The suggested 
beaming (OPERA like experiment) of neutrinos from CERN (or FNL) across the Earth to 
IceCube might provide a novel signature for muon and anti-muon neutrino mixing while better 
focusing flavor parameters . Nevertheless the present atmospheric neutrino noise that should 
arise in Deep Core data is also bearing important records on neutrino mixing and eventual CPT 
violation. Indeed we used previous (Deep Core) expected rate to foresee (over the same expected 
spectra) deviation due to CPT (MINOS like) violations. Nevertheless a more accurate study 
did lead us to foresee a different spectra [8J , to be elaborated here, more in detail, based on the 
known Super-Kamiokande detection, scaled to Deep Core volumes (and environment). Earlier 
[8j preliminary result is here more deepened and conflrmed, starting from primordial atmospheric 
spectra, source of atmospheric neutrinos up to their interactions and muon event rate whose 
tracks are projected along each channel string. Results are somehow surprising respect earlier 
ones. 



2. Prom the atmospheric cosmic ray to the event rate in Deep Core 

Let us remind the neutrino energy flux deflned as 
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that we parameterized, here for the flrst time, fltting the well known curve from refs.[2l [3] 
for muon neutrinos with the following analytical formula: 
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where the energy is expressed in GeV unity; the energy range of this formula embrace 
the hundred MeV up to TeV. This averaged isotropic spectra must be enhanced (empirically) 
at horizons (slightly at low, GeVs energy and more at high TeV energy) because of longer 
pion, kaons and muons decay flight. We disregarded for Deep Core and IceCube at South 
Pole the geomagnetic cut off influence; (its role at lowest energy has been taken into account 
at SuperKamiokande for azimuthal directions). To take into account of the enhancement (of 
horizontal respect to vertical rate) we derived (to flt the SK data) a correction factor to be 
multiplied to previous flux: 
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where the energy is in GeV unity, and horizontal angle (complemental to zenith angle z) 
range from (horizon) to 7r/2 (upward vertical direction). The total mass target in Deep Core 
for incoming neutrino is taken from the effective volume estimated (for each energy) in [4J, 
without online veto. The ratio between all (fully contained (FC) with partially contained (PC) 
and upward stopping and through going, i.e. without on line veto) events versus only (FC+PC, 
i.e. on-line veto) is unity up to ^ 100 GeV, reaching a ratio ^ 1.5 at 400 GeV, or ~ 2 at 
TeV. In these highest energies the channel rate is low and it is not so relevant in our neutrino 
oscillation study, mostly focused below 200 GeV and below maximal 50 channels. Therefore we 
considered here all the events without veto, almost coincident with the ones with veto within 
50 channels count, within the following Deep Core energy dependent (as expected following ^ 
in GeV energy unity), threshold volume: 
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Neutrino Charge Current (CC) cross section for muon production is considered (see ref. [^) 
to be linear increasing with energy, in the energy range considered, while antineutrino cross 
section is taken about a half of the neutrino one. Electron neutrino and tau neutrinos CC and 
all flavor neutral current (NC) event are source of spherical showers whose role may influence 
(and pollute) only a few channel (3-5) and are only mentioned here: they will be discussed in 
detail elsewhere. 

The calculation includes neutrino oscillation probability as a function both of energy and 
arrival observation angle, given the earth chord distance for each angle below the horizon, 
Piy^^iy.(Eiy,0),i = 1,2,3; it means that we take into account three exact neutrino flavor 
oscillation (see [7J), discussing the out-coming (observable) muon tracks. We took into account 
also the matter influence on the mixing inside the Earth assuming at each zenith angle the 
corresponding average Earth density. We will show elsewhere that the exact integral in variable 



Earth density do not change our present results at 10-100 GeV ranges, but mostly induce 
changes at lower energies. 

We obtain thus the differential event number (per unit energy and solid angle and unit time) 
in function of neutrino energy and the horizontal angle 9, in our definition starting from 0° at 
horizon, equal to 90° at vertical upward direction: 
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Event rate in one year of observation is obtained integrating over 9 from 37° below the horizon 
(equal to zenith angle cos(z) = cos(7r/2 — 9) ^ sin(6>) = 0.6 in Deep Core zenith angle) to 
90°, and integrating over an appropriate energy bin: 

-^<-) - =^<-' (^) - /r^7:::^<--' ■ m 

where Na is the Avogadro number (multiplied by 10^, because ice density picE ^ 10^ kg m~^) 
and V{Ejj) is the effective interaction volume in Deep Core (expressed in m^, like eq. (jl])). 
Resulting events distribution are shown in figs. El [5] and El for ideal narrow and real wide energy 
width. Later on the same spectra of muon tracks has to be projected along the vertical axis 
of Deep Core, as shown finally in fig. [71 The muon born in Charged Current interaction has 
nearly half the neutrino energy, and undergoes an average energy loss of dEj^/dX 2^ 200 MeV 
m~^ or 1/5 GeV m~^. Therefore the muon length (keeping in mind that kinematically holds 
Ejj, ^ Ey/2) in water is approximately: — bEy/2 m, with energy as usual expressed in 
GeV. We used indeed the exact linear-logarithmic growth in water: 

= 3.9-1 . (1 + 1^ . 10-3 E,)j m. (7) 

To obtain a comparable plot of events distribution versus channel number we assumed an average 
projection factor, (cos(7r/2 — 9)) = (sin6>) of muon track along the DOM strings due to the 
spread of the muon track direction (in range 0.6 < sin0 < 1) first to be 0.8 at same plane; 
however an additional projection factor 0.8 is needed because muons tracks may be often skew 
respect the string hne. The overall factor ranges (in average between 0.8^ and 0.8). We assumed 
the medium projection factor equal to 0.72. 
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Each channel to channel distance corresponds no longer to 2.8 GeV (as for a vertical muon) 
but because of the inclined muon track projected along the vertical string < sin9 0.72, 
on average, to a wider energy band of 2.8/0.72 = 3.88 GeV. This projecting factor smear the 
sharp decay on rate event due to the oscillation and it mask most of the low energy neutrino 
oscillation, whose inner traces has been revealed by SK only in last decade. However the most 
severe deepening in the spectra (see our result versus earlier one in fig. [7j) at about 20 GeV 
survive the smearing and it marks our expected count rate. 



3. Conclusions 

The Deep Core telescope rate should be soon pubhshed: their spectra (projected and counted 
along the channel) might be foreseen now and compared soon with the records. Even if we did 



first used the Deep Core predictions [7J, in late papers [8J we did wonder and did suggest a 
different expected rate. In the present paper we show step by step the procedure to foresee the 
Deep Core rate based on known Cosmic Ray rate and on SK muon neutrino flux, flux interacting 
in Deep Core (variable) volume, leading, after mixing (see Appendix) to survived muon (and 
anti-muons) neutrinos whose lepton tracks may rise in Deep Core Optical Module string. Our 
results differ from other ones ^ , see fig [HI we tried to verify the estimates in different check 
(versus other experiments) and we believe that they are (at first approximation) correct. The 
Deep Core rate shape, its intensity might be soon verified and tested. If the expected (CPT 
conserved estimate ) are met, Minos claim might be rejected. If at some channels (mostly 3- 
4-5) there will be an over production of signals, than the Minos CPT violated hint may be 
taken seriously. It should be taken into account that NC (neutral current) events at tens GeV 
may spread light in spherical shape increasing and polluting the (3-6) channel rate; these a few 
thousands (2^ 3000 — 4000) events may enrich channel 2-5; however their timing structure and 
their uneven light intensity ratio may probe their different nature (see [8]) from the up-going 
vertical muon ones, possibly to be separated easily. Therefore these new Deep Core records may 
probe at minimum channel gate (6-8) an atmospheric suppression offering a rare windows for 
20 GeVs neutrino astronomy. To observe such an astronomy in the whole arrival direction (not 
just along the zenith but also in its azimuth) we suggest the doubling of the present Deep Core 
array strings, allowing a contemporaneous double string detection. The discovery of twenty GeV 
neutrino astronomy in Deep Core may find a bonus also in rare, but nearby and powerful. Solar 
Flare neutrino brightening [10\ 

4. Appendix: 3 flavor neutrino mixing 

The evolution equation for three neutrino mixing in matter reads : 

i^^a = H^a, with (UMU^ + ^) , a = 1, 2, 3 (9) 

which is a Schrodinger equation with a matter potential A, 
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where Acc = 2^/2EGFNe = 0.76 • 10"^el/2(^)(-^) is the potential only due to Ue 
charged current interaction. 

This equation is a set of 9 differential equations whose solutions for probability conversion 
Pa(3 — with a,/3 = e,iJ.,r, has boundaries conditions Pa(3{0) equal to 1 for the starting 

flavor /X, for the others (electron conversion into muon is negligible). This system of equations 
is symmetric since P^^ = P/Sa- 

Among the simplest numerical solutions to show, we found a compact survival probability, for 
instance in vacuum for muon to muon, at a distance equal to the Earth diameter, corresponding 
to = 90°, as a function of neutrino energy: 

P^^ = 1 0.240128 + 0.304714 e-2-5859i/£;n ^ 0.455158 e-^^-^^oW^^I^ 

The more general cases taken into our study, with earth matter density included, are 
expression too wide and too complex (few pages of mathematica program) to be written here. 
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Figure 1. The differential number event count a year in Deep Core as a function of the muon 
neutrino energy and angle, observable in narrow scale (top, 3-35 GeV) assuming the vacuum 
among the neutrino flight). The horizontal axis is the energy in GeV, the vertical one is the event 
rate a year in Deep Core, the depth axis is the zenith angle variable in radiant. 
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Figure 2. As above the differential number event count a year in Deep Core as a function 
of the muon neutrino energy and angle, observable in wide scale (top, 3-65 GeV) assuming an 
average Earth density p — 1 during up-going neutrino flight within Earth. The horizontal axis is 
the energy in GeV, the vertical one is the event rate a year in Deep Core, the depth axis is the 
zenith angle variable in radiant. Most of the inner oscillating structure is lost because of the 
Deep Core un-ability to reveal a few GeVs muon track. The last longest flux suppression at 20 
GeV survive and is source of a spectra deepening in Channel 6-9 
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Figure 3. The ideal count rate 
that Deep Core would be able to 
reveal on the muon track length if 
it could disentangle muons within 
20 MeV energy width. It is clear 
the marked difference between non 
oscillated case (red non oscillat- 
ing curve), the mixed CPT con- 
served case ( oscillating blue curve) , 
and the thin-dashed CPT violated 
(green) alternative signal. This is 
in vacuum case. 
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Figure 4. As previous figure, 
where neutrino mixing occurs in 
matter. Note the little difference in 
spectra at the very low energies. 
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Figure 5. As previous figure the averaged 
(within 1 GeV) spectra of atmospheric 
neutrino if Deep Core would be able to 
disentangle within 1 GeV energy width. As 
in previous curves, the non oscillated one 
is red and non oscillating , the mixed CPT 
conserved case is still oscillating and it is 
a blue curve, while the thin-dashed CPT 
violated (green) curve is an alternative 
signal. Most, but not all, of the low energy 
oscillations are lost. 



Figure 6. A possible real count rate as a 
function of the energy within a 2.8 GeV 
width (as each channel from a signal of 
vertical muon in Deep Core). As above 
the red curve marks the un-oscillated case, 
the blue one the CPT conserved version 
and the thin-dashed green line stands for 
the CPT violated solution. The smearing 
of the earliest oscillation is complete. Few 
mixing structure remains at low channels. 
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Figure 7. The most probable realistic rate count rate a year as a function of the channel within 
a 3.88 GeV width projected within a factor 0.72 (as the channel distance in Deep Core). As 
above the red curve marks the un-oscillated case, the blue one the CPT conserved version and 
the thin-dashed green line stand for the CPT violated solution. The smearing of the earliest 
oscillation is hidden. But the deepening structure remains at 6-9 channels remind of the muon 
disappearance at twenty GeVs. The depression of the muon flux in this region is nearly a factor 
ten, making this window silent and ideal to search for a Neutrino Astronomy 
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Figure 8. Our expected rate count rate (histograms) for upgoing muon tracks versus the most 
recent Deep Core expectations (dots) As above the red curve marks the un-oscillated case, 
the blue one the CPT conserved version and the thin-dashed green line stand for the CPT 
violated solution. The discrepancy is relevant and we cannot find an explanation for it. Only 
above channel 50 (where we cannot await for measures) the two rate maybe comparable; in all 
the channel range our estimate differ greatly in value and shape from the [4J. The presence 
of additional neutral current shower events will increase channel 2-5 (almost doubling their 
number), but they maybe in principle disentangled from muon upward tracks by their uneven 
intensity ratio and their timing signature [S]. 



